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Preface 

The  purpose  of  this  study  was  to  develop  a  computer  model  of  the  dynamics  and 
controls  for  a  robot  spacecraft  used  to  capture  a  tumbling  and  spinning  satellite.  The 
robot  spacecraft  is  to  detumble  and  despin  the  other  spacecraft  so  it  can  be  repaired  and 
returned  to  useful  service.  I  would  like  to  acknowledge  the  help  of  Lt  Col  Joseph  W. 
Widhalm  in  getting  me  started  on  this  endeavor,  which  was  a  follow-on  to  his  previous 
work.  I  would  also  like  to  acknowledge  the  help  of  Dr.  Curtis  H.  Spenny  who  later  filled 
in  as  my  thesis  advisor  after  Lt  Col  Widhalm’s  departure  in  July  of  this  year. 
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unit  vector  set  arbitrarily  fixed  at  the  center  of  mass  of  Body  0  (OMV) 
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nongravitational  external  force  on  body 

interaction  force  on  body  a  transmitted  through  joint  j 
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i,] ,fc  =  unit  vector  set  representing  inertial  planetocentric  coordinate  system 
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unit  vector  set  arbitrarily  fixed  at  the  center  of  mass  of  the  target 
satellite  and  aligned  with  its  principal  axes 

number  of  rotational  degrees  of  freedom  in  the  n-body  system 
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vector  from  the  system  center  of  mass  to  the  center  of  mass  of  body  a 

nongravitational  external  torque  on  body  a 

constraint  torque  on  body  a  at  joint  j 

spring  damper  torque  on  body  a  at  joint  j 

time 
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nine-element  control  vector 
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=  nine-element  state  vector 
=  augmented  thirteen-element  state  vector 
=  three-element  joint  translational  vector 
=  unit  dyadic 

=  planet’s  gravitational  constant 
=  angle  of  rotation  about  axis 
=  coning  angle  of  the  target  satellite 

=  planetocentric  position  vector  of  n-body  system  center  of  mass 
=  unit  vector  in  the  direction  of  p 

=  planetocentric  position  vector  to  the  center  of  mass  of  body  x 
=  inertia  dyadic  of  body  x 
=  pure  spin  rate  of  target  satellite 
=  precession  rate  of  the  target  satellite 
=  angular  velocity  of  body  0 
=  angular  velocity  of  body  x 
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Abstract 

The  problem  of  detumbling  and  despinning  a  freely  spinning  and  processing 
axisymmetric  target  satellite  using  an  orbital  maneuvering  vehicle  is  considered.  The 
axisymmetric  orbital  maneuvering  vehicle  is  equipped  with  a  multibody  grappling  arm 
assembly  to  capture  the  target.  Counter- masses  are  used  to  maintain  dynamic  balancing 
and  stability  throughout  the  deployment  of  the  arm  and  the  subsequent  spin-up  of  the 
grappling  device  prior  to  docking.  The  five-body  system  is  modeled  using  Eulerian- based 
equations  of  motion  developed  by  Hooker  and  Margulies.  Open- loop  control  laws  are 
formulated  to  deploy  the  grappling  arm  assembly  and  spin-up  the  grappling  device  using 
internal  motor  torques.  A  Liapunov  technique  is  applied  to  derive  a  nonlinear  feedback 
control  law  that  drives  the  docked  system  to  a  final  spin-stabilized  state  of  equilibrium. 
External  thrusters  are  used  to  maintain  the  absolute  motion  of  the  system  during  this 
process.  Variations  in  grappling  arm  length,  target  coning  angle,  and  response  times  are 
examined  for  design  purposes.  State  and  control  histories  are  presented  and  the  results 
from  this  five-body  model  are  compared  with  the  Widhalm  and  Conway  two-body  study. 
The  simulation  indicates  that  the  required  control  magnitudes  are  higher  for  the  five-body 
model  but  are  still  quite  reasonable.  The  addition  of  the  grappling  arm  assembly  adds 
both  realism  and  flexibility  to  the  capture  problem. 


REMOTE  ORBITAL  CAPTURE  USING  AN  ORBITAL  MANEUVERING  VEHICLE 
EQUIPPED  WITH  A  MULTIBODY  GRAPPLING  ARM  ASSEMBLY 


I.  Introduction 

The  servicing  and  repair  of  satellites  in  orbits  beyond  the  direct  reach  of  the  Space 
Shuttle  requires  a  teleoperator  spacecraft  or  orbital  maneuvering  vehicle  (OMV)  to  dock 
with  these  target  satellites  and,  if  necessary,  return  them  to  an  accessible  orbit.  If  the 
target  has  experienced  a  control  system  failure,  it  may  be  necessary  to  detumble  it.  If  the 
target  satellite  is  spin  stabilized  it  may  be  necessary  to  despin  it.  Docking  followed  by 
despinning  or  detumbling  is  defined  here  as  remote  orbital  capture. 

Previous  work  by  Widhalm  and  Conway  (9)  considered  both  the  detumbling  and 
despinning  aspects  of  remote  orbital  capture,  and  the  requirement  to  control  the  absolute 
motion  of  the  two-body  system  using  feedback  control.  In  their  model  (Fig.  1),  an 
axisymmetric  OMV  was  docked  with  a  freely  spinning  and  processing  axisymmetric  target. 
The  OMV  was  equipped  with  a  two  rotational  degree  of  freedom  joint  that  connected 
with  a  fixed  appendage  on  the  target  body.  This  joint  was  free  to  translate  along  the 
surface  of  the  OMV  to  provide  positioning  adjustments  during  docking.  Despinning  and 
detumbling  were  accomplished  using  feedback  control  by  driving  the  joint  to  the  OMV’s 
axis  of  symmetry  and  applying  internal  motor  torques  to  the  target  through  the  connecting 
joint.  While  this  was  being  done,  the  OMV  external  thrusters  were  fired  to  control  the 
absolute  motion  of  the  two-body  system.  It  was  shown  that  this  "detumbling  process  is 
quite  benign  and  that  the  required  control  magnitudes  are  small”  (9:657). 

One  key  assumption  made  in  the  Widhalm  and  Conway  model  was  that  the  target  had 
a  fixed  appendage  that  the  OMV’s  joint  connected  with  to  form  the  two- body  docked 
configuration.  In  reality,  however,  a  target  satellite  would  probably  not  have  this 
convenient  appendage  to  serve  as  a  docking  point.  What  would  really  be  required  is  for 
the  OMV  to  be  equipped  with  a  grappling  arm  assembly  which  it  could  deploy  to  capture 
and  detumble  the  target  satellite. 

This  thesis  extends  the  Widhalm  and  Conway  two-body  model  by  taking  into  account 
the  requirement  for  a  grappling  arm  assembly  attached  to  the  OMV.  The  new  model 
consists  of  five  constant  mass  rigid  bodies.  Simulation  of  the  capture  process  now 
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becomes  a  three  step  rather  than  a  single  step  process.  The  first  step  is  the  deployment  of 
the  grappling  arm  assembly  up  to  but  not  including  docking  with  the  target  vehicle.  The 
deployment  of  the  grappling  arm  significantly  alters  the  dynamics  of  the  system,  and 
one  important  aspect  of  the  problem  is  to  maintain  dynamic  balancing  and  stability 
throughout  the  process..  The  second  step  is  the  spin-up  of  the  grappling  device  to  match 
the  pure  spin  rate  of  the  target  ,  just  prior  to  docking.  The  actual  docking  process  itself 
is  not  treated  in  this  thesis,  and  is  handled  as  a  discrete  event.  The  final  step  once  the 
target  satellite  is  joined  with  the  OMV  is  to  apply  the  required  internal  motor  torques 
and  OMV  thruster  firings  using  feedback  control  to  detumble  and  despin  the  target  and 
drive  the  entire  system  to  a  spin-stabilized  state  of  equilibrium. 

Both  the  system  initial  configuration  and  the  desired  final  state  are  defined.  The 
applicable  equations  of  motion  for  the  five  body  system  are  then  derived,  and  the 
required  control  laws  are  formulated  for  each  of  the  three  steps  in  the  capture  process. 
The  results  of  the  simulation  show  the  magnitude  of  both  the  internal  and  external  control 
torques  required  to  achieve  the  desired  outcome.  This  information  would  be  necessary  in 
eventual  hardware  design  for  sizing  motors  and  thrusters,  and  in  performing  structural 
analysis  of  the  system. 
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System  Configuration 


The  OMV  model  for  this  thesis  consists  of  five  rigid  bodies  as  depicted  in  Fig.  2. 
These  five  bodies  are;  the  original  axisymmetric  reference  or  base  body  used  in  the 
Widhalm  and  Conway  model,  the  grappling  arm,  the  counter-mass,  and  the  grappling 
device  with  its  corresponding  counter-mass.  The  ei,e2,e3  coordinate  system  used  here  is 
a  rotating  system  fixed  at  the  center  of  mass  of  Body  0,  not  at  the  center  of  mass  of  the 
system.  The  e2  axis  is  aligned  along  the  direction  of  the  translating  joint  and  the  the 
grappling  arm  assembly,  and  the  ea  axis  is  the  axis  of  symmetry  for  the  reference  body. 
There  is  no  specified  coordinate  system  fixed  at  the  center  of  mass  of  the  system  and  all 
equations  are  eventually  expressed  in  the  e],,e2,e3  system. 

Body  0,  the  reference  or  base  body,  has  angular  velocity  components  wqi,  wo2,  wqs 
corresponding  to  its  three  principal  axes,  ei,e2,e3.  Also,  Body  0  contains  the  three 
external  variable  control  thrusters.  The  three  external  control  torques  generated  ,  Tl,  T2, 
T3,  again  correspond  to  the  ei,e2,e3  principal  axes.  To  simplify  the  problem,  it  is 
assumed  that  firing  the  thrusters  does  not  affect  the  mass  properties  of  the  OMV  and  the 
entire  model  is  assumed  to  be  a  constant  mass  system. 

There  are  six  additional  rotational  degrees  of  freedom  for  the  overall  OMV  system, 
specified  by  their  respective  axes  of  rotation,  gi,  g2,  gs,  g^,  gs,  ge,  for  a  total  of  nine 
rotational  degrees  of  freedom.  The  amount  of  rotation  about  each  axis  is  specified  by 

Y\-  yj-  Ys.  Y*-  Ys .  Yt  ,  respectively.  There  are  also  six  internal  motor  torques,  TGI, 
TG2,  TG3,  TG4,  TG5,  TG6  ,  again  corresponding  to  the  six  axes  of  rotation. 

Body  1  is  connected  to  the  reference  body  by  Joint  1  which  is  the  same  type  of 
translating  joint  used  in  the  Widhalm  and  Conway  model.  Joint  1  has  two  rotational 
degrees  of  freedom,  one  about  gj  which  is  parallel  to  the  ei  axis,  and  one  about  g2  which 
is  aligned  with  the  grappling  arm.  The  amount  of  rotation  of  Body  1  relative  to  Body  0 

in  the  positive  ei  direction  is  measured  by  the  angle  Yi.  Body  1  is  the  grappling  arm 
itself,  and  does  not  include  the  two  end  masses.  Body  3  and  Body  4,  which  represent  the 
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Body  2 


Joint  1  Body  1 


BodyO 
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grappling  device  assembly.  However,  Body  1  does  include  the  mass  of  the  translating 
joint  assembly.  Joint  1. 

Body  2  is  the  counter-mass  which  is  designed  to  balance  out  the  combined  effects  of 
Body  1,  Body  3,  and  Body  4  so  that  no  cross  products  of  inertia  terms  are  generated 
during  the  deployment  phase.  Body  2  is  connected  to  Body  1  with  Joint  2.  Its  only 
degree  of  freedom  is  about  the  ga  axis  which  is  parallel  to  the  ei  axis.  The  amount  of  this 

rotation  relative  to  Body  1  is  measured  by  and  is  in  the  negative  ei  direction. 

Body  3  is  the  grappling  device.  Note  that  it  is  not  the  scope  of  this  thesis  to  design 
or  specify  an  actual  hardware  device.  Body  3  is  merely  a  mathematical  representation  of 
the  grappling  device.  Similarly,  the  actual  docking  process  is  not  part  of  this  simulation, 
and  is  merely  a  discrete  event.  The  system  is  still  undocked  at  the  end  of  the  deployment 
and  spin-up  phases  and  already  docked  at  the  start  of  the  detumble/despin  phase. 

Body  4  is  another  counter-mass  with  exactly  the  same  mass  and  inertia  properties  as 
Body  3.  Its  purpose  is  to  spin  at  an  equal  but  opposite  rate  as  Body  3  during  the  spin-up 
phase.  This  will  cancel  out  the  moment  created  by  spinning  up  Body  3. 

Body  3  and  Body  4  are  both  attached  to  Body  1  by  Joint  3  which  must  allow  for  three 
rotational  degrees  of  freedom.  It  allows  Body  3  and  Body  4  to  rotate  as  one  body  about 
the  g^  axis  which  is  parallel  to  the  ei  axis.  The  amount  of  this  rotation  is  measured  by 

the  angle  and  the  direction  is  in  the  negative  ei  direction.  Note  that  >4  is  defined  to  be 
in  the  0°  position  when  the  g5  and  ge  axes  are  perpendicular  to  Body  1 ,  as  shown  in  Fig.  2. 
In  addition  there  is  the  spin  rotation  of  Body  3  about  the  gg  axis,  and  the  spin  rotation  of 
Body  4  about  the  gg  axis.  Note  that  gg  and  gg  are  the  same  axis,  with  the  spin  of  Body  3 
in  the  positive  direction  and  the  spin  of  Body  4  in  the  negative  direction.  Note  also  that 
Body  3  and  Body  4  have  axisymmetric  inertia  properties  with  their  axis  of  symmetry 
being  this  gg/gg  axis. 

In  the  initial  configuration  prior  to  deployment  of  the  arm,  the  OMV  is  positioned 
relative  to  the  target  satellite  such  that  the  OMV’s  axis  of  symmetry,  eg,  is  parallel  to  the 
target’s  angular  momentum  vector.  The  target  satellite  is  both  tumbling  and  spinning, 


where 


y  T  ’  V  T  ■>  ^  T 
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correspond  to  the  target’s  coning  angle,  precession  rate,  and  spin  rate  (about  its  ns  axis), 
respectively.  These  three  properties  are  related  as  shown  by  Greenwood  (5:386)  and 
repeated  here: 


V 


ii^^cosy^. 


where  Ixx  is  the  target’s  moment  of  inertia  about  the  ni  and  ns  axes,  and  is  the  target’s 
moment  of  inertia  about  its  axis  of  symmetry,  ns  .  The  OMV  is  positioned  such  that  the 
target’s  center  of  mass  is  lined  up  on  the  OMV’s  es  axis.  The  OMV  is  spun- up  so  that  its 
angular  velocity  component  wqs  about  the  63  axis  exactly  matches  the  target’s  precession 
rate,  while  the  wqi  and  wqs  components  remain  zero.  The  OMV  is  therefore  in  a  state  of 
pure  spin.  Besides  the  target’s  pure  spin  rate,  the  OMV  and  the  target  will  have  zero 
relative  motion.  The  OMV  is  also  positioned  at  a  pre-determined  stand-off  distance  from 
the  target,  such  that  when  Joint  1  translates  a  specified  amount  in  the  63  direction  and  the 

grappling  arm  is  raised  by  a  corresponding  angle  y ,  ,  the  grappling  device  will  be  able  to 
successfully  dock  with  the  target  as  shown  in  Fig.  3.  This  stand-off  distance  depends  on 
the  target’s  length  and  coning  angle,  the  length  of  the  grappling  arm,  and  the  allowable 
amount  of  joint  translation.  It  is  assumed  that  the  complete  initial  state  of  the  target  has 
been  accurately  determined  by  the  OMV’s  computer,  that  the  OMV  has  been  properly 
placed  into  the  required  position,  and  has  already  been  spun-up  to  match  the  target’s 
precession  rate.  It  is  with  this  initial  configuration  that  this  thesis  begins. 

The  first  phase  of  the  simulation  is  the  deployment  of  the  grappling  arm  assembly.  In 
the  initial  state  Bodies  1,  2,  3,and  4  are  stored  horizontally  on  the  surface  of  the  OMV  so 

that  Yi .  Va.  74  are  all  zero.  Joint  1  is  commanded  to  begin  translating  in  the  63 
direction  toward  the  es  axis  of  symmetry,  and  internal  motor  torques  are  applied  to  begin 
erecting  the  grappling  arm  assembly.  It  is  the  intent  of  this  thesis  to  perform  the  entire 
deployment  process  and  still  maintain  the  OMV  in  a  state  of  pure  spin  at  the  constant 
initial  rate  with  a  minimum  application  of  external  control  torques,  if  any.  It  is  desired 
not  to  fire  these  external  thrusters  unless  absolutely  necessary  in  order  to  conserve  the 
attitude  maneuvering  fuel. 
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It  is  therefore  critical  in  order  to  maintain  dynamic  balancing  and  stability,  that;  1) 
the  combined  center  of  mass  of  Bodies  1,  2,  3,  and  4  always  lies  on  the  63  axis,  and  2)  the 
inertia  properties  of  the  counter-mass  Body  2  are  such  that,  when  combined  with  the 
inertia  properties  of  Bodies  1,  3,  and  4,  the  ea  axis  is  always  the  axis  of  symmetry  for  all 
five  bodies  so  no  cross  products  of  inertia  terms  are  generated.  To  solve  this  problem,  the 
mass  properties  of  Body  1,  Body  3,  and  Body  4  must  be  such  that  the  total  mass  center  of 
these  three  masses  is  at  Joint  2  ,  which  is  located  at  the  geometric  center  of  the  grappling 
arm  and  is  also  the  center  of  mass  of  Body  2.  Therefore,  the  combined  center  of  mass  of 
Bodies  1,  2,  3,  and  4  is  located  at  Joint  2  and  the  control  problem  now  is  to  maintain 
Joint  2  always  along  the  63  axis  during  the  deployment  process.  This  is  achieved  with  the 
use  of  an  open-loop  control  algorithm  and  the  application  of  internal  motor  torques.  It  is 
also  desired  to  lock  the  rotations  about  g2,  g4,  gs,  and  ge  fixed  during  this  process,  since 
they  are  not  required.  As  for  the  cross  product  of  inertia  problem,  Body  2  should 

maintain  a  constant  relationship  with  Body  1  such  that  Vs  is  always  exactly  twice  as  much 

as  Vi,  but  in  the  opposite  direction.  Again,  this  is  handled  with  the  appropriate  open-loop 
algorithm  and  the  application  of  the  required  internal  motor  torques. 

Throughout  this  process  if  the  thrusters  are  not  fired  the  OMV’s  system  mass  center 
will  not  move  in  inertial  space  since  no  external  forces  or  torques  are  being  applied. 
However,  the  system  mass  center  will  be  moving  "upward"  relative  to  the  OMV’s  base  and 
so  the  center  of  mass  of  the  OMV’s  base,  the  reference  Body  0,  will  translate  "downward" 
in  inertial  space  along  the  vertical  es  axis.  The  reference  body’s  position  relative  to  the 
target  will  therefore  be  changing  during  the  process.  Although  this  translation  can  easily 
be  calculated,  it  will  be  assumed  for  the  purposes  of  this  thesis  that  the  OMV  has  a 
separate  control  system  that  will  compensate  for  it  by  applying  the  appropriate 
translational  thrusters,  thereby  maintaining  zero  relative  motion  between  the  base  of  the 
OMV  and  the  target. 

The  arm  will  be  deployed  up  to  but  not  including  docking.  Body  3  will  then  be 
assumed  to  be  at  an  infinitely  small  distance  away  from  the  target.  In  reality,  this  would 
probably  be  the  best  opportunity  to  measure  the  spin  rate  of  the  target  accurately  using  a 
sensor  on  Body  3.  At  this  point  the  second  phase  of  the  simulation  begins.  Body  3  will 
be  spun- up,  from  rest,  to  the  same  rate  as  the  pre-determined  spin  rate  of  the  target.  At 
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the  same  time,  Body  4  will  be  spun-up  at  an  equal  and  opposite  rate  ,  cancelling  out  the 
moment  created  by  Body  3,  in  order  to  maintain  dynamic  balance  and  stability.  Body  3 
will  now  have  zero  relative  velocity  in  relation  to  the  target,  and  will  be  ready  to  grab 
hold  of  it.  Prior  to  docking,  the  entire  OMV  system  should  be  in  a  state  of  pure  spin 
about  the  es  axis  without  the  application  of  the  external  thrusters. 

The  third  and  final  phase  of  the  simulation  begins  immediately  after  docking  with  the 
target.  After  this  discrete  event,  the  total  system  is  still  treated  as  a  five-body  problem 
by  simply  adding  the  mass  and  inertia  properties  of  the  target  to  Body  3,  and  treating  the 
target  and  grappling  device  as  one  single  body  (still  labeled  Body  3  -  See  Figs.  2  and  3). 
The  closed-loop  feedback  system  is  now  turned  on,  and  both  internal  motor  torques  and 
external  thruster  torques  are  applied  to  drive  the  system  to  the  final  desired  state.  This 

final  state  is  that  in  which  Joint  1  is  centered  at  the  63  axis,  Vi  =  90°  ,  Vs  =  180°  ,  and 

>4  =  90°  so  that  Body  1,  Body  2,  Body  3  (including  the  target),  and  Body  4  are  all  in  the 
vertical  position,  as  shown  in  Fig.  4  and  the  entire  system  is  in  a  spin-stabilized  state  of 
equilibrium. 
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Body  2  Requirements 


As  discussed  in  the  previous  section.  Body  2  with  its  center  of  mass  at  Joint  2  must 
have  the  equivalent  inertia  properties  as  Body  1,  Body  3,  and  Body  4  with  their  combined 
mass  center  also  at  Joint  2.  Body  2’s  inertia  requirements  can  be  determined  by  a  simple 


application  of  the  parallel  axis  theorem. 
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(i),  (i)  refer  to  the  principal 

mass  moments  of  inertia  of  each  body. 

Note  that  when  all  bodies  are  initially  in  the  stored  position,  these  principal  moments  of 
inertia  are  parallel  to  the  ei,  ej,  es  axes,  respectively.  Note  that  L12  is  the  length  from 
the  center  of  mass  of  Body  1  to  Joint  2,  and  L13  ,  L14  are  both  the  length  from  the  center 
of  mass  of  Body  1  to  Joint  3.  Also,  L31  is  the  length  from  the  center  of  mass  of  Body  3 
to  Joint  3  which,  prior  to  docking,  is  equal  to  L41  . 


w. 
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The  Hooker-Margulies  Equations 

The  appropriate  equations  of  motion  now  need  to  be  developed  to  describe  the 
five-body  model.  Fletcher,  Rongved,  and  Yu  (4)  derived  the  dynamical  attitude  equations 
for  a  two-body  satellite.  Hooker  and  Margulies  (7)  then  generalized  these  equations  for 
an  n-body  satellite.  They  used  Newton’s  equations  to  eliminate  the  moments  of  the 
unknown  interaction  forces  (i.e.  the  moments  of  the  reactive  forces  transmitted  through 
the  joints)  from  the  Euler  equations.  They  therefore  derived  a  complete  set  of  3n  scalar 
equations  for  an  n-body  system  that  are  free  of  unknown  joint  constraint  forces  but  still 
contain  the  unknown  constraint  torques.  There  are,  however,  two  restrictions  placed  on 
these  equations.  First  of  ail,  the  topology  of  the  overall  configuration  of  bodies  and 
interconnecting  joints  must  be  equivalent  to  a  topological  tree  (i.e.  no  closed  loops). 
Second,  the  relative  motion  about  a  joint  is  assumed  to  be  rotational  and  does  not  allow 
translation  of  either  body  relative  to  the  joint.  (7:123) 

Hooker  (6)  in  a  follow-on  paper  showed  that  the  the  constraint  torques  could  also  be 
eliminated.  This  reduces  the  number  of  dynamical  attitude  equations  from  3n  to  r  where 
r  is  the  number  of  rotation  degrees  of  freedom  for  the  system.  The  number  of  dependent 
variables  is  also  reduced  from  the  set  of  3n  angular  velocity  components  to  3  angular 
velocity  components  of  the  reference  body  and  an  additional  r-3  relative  angular 
rotational  rates.  The  original  Hooker-Margulies  equations  are  written  for  all  the  bodies 
lying  on  one  side  of  a  selected  joint  and  subsequently  added.  The  interaction  torques  all 
cancel  in  pairs  with  the  exception  of  the  constraint  torque  at  the  selected  joint.  If  this 
particular  joint  has  a  rotational  degree  of  freedom  about  an  axis  g,  then  the  dot  product 
of  g  and  the  expression  for  the  constraint  torque  is  zero,  and  so  writing  the  dot  product 
and  setting  it  to  zero  yields  an  equation  that  is  free  of  the  constraint  torque.  This  process 
can  then  be  repeated  for  each  degree  of  freedom  at  each  joint  to  eliminate  all  the 
unknown  constraint  torques  from  the  equations.  (6:1205-1207)  This  subsequent  set  of  r 
equations  for  an  n-body  system  with  r  rotational  degrees  of  freedom  is  called  the 
modified  Hooker-Marqulies  equations. 

An  equivalent  set  of  dynamical  attitude  equations  of  motion  could  be  derived  using 
the  Lagrangian  method.  The  Lagrange  approach  to  deriving  the  equations  of  motion  has 
the  advantage  that  the  constraint  torques  never  appear,  and  that  the  number  of  equations 
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is  r,  the  number  of  rotational  degrees  of  freedom  for  the  system.  However,  the  resulting 
equations  would  not  be  written  in  terms  of  physical  body  axes,  as  are  the  HM  and 
modified  HM  equations.  As  a  result,  it  would  be  exceeding  difficult  to  modify  the 
Lagrange  equations  in  order  to  adapt  them  to  an  active  control  system  and  to  include  the 
the  effects  of  joint  motion.  (6.1205) 

One  other  important  advantage  of  the  modified  HM  equations  is  that  they 
significantly  reduce  the  amount  of  computer  time  by  only  having  to  integrate  numerically 
r  equations  rather  than  3n  equations.  (6:1205) 

As  previously  mentioned,  neither  the  HM  equations  nor  the  modified  HM  equations 
account  for  translation  of  the  joints  relative  to  the  bodies  adjacent  to  the  joint.  In  this 
thesis,  as  in  the  Widhalm  and  Conway  model.  Joint  1  connecting  Body  0  and  Body  1 
translates  relative  to  Body  0,  the  reference  body.  This  was  dealt  with  by  Conway  and 
Widhalm  (1)  with  an  extension  of  the  modified  HM  equations  that  permits  the  translation 
of  the  joint.  It  is  therefore  this  extended  version  of  the  modified  HM  equations  that  can 
now  be  applied  to  this  thesis. 

Even  though  the  modified  HM  equations  eliminate  both  the  constraint  forces  and  the 
constraint  torques,  this  information  may  later  be  needed  in  performing  a  structural 
analysis  of  the  system.  Fortunately,  it  is  relatively  easy  to  recover  this  information  as 
shown  in  references  (6)  and  (7). 
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Application  of  the  Hooker-Margulies  Equations 

The  dynamical  attitude  equations  for  the  5-Body  OMV  system  (n=5)  with  its  9 
rotational  degrees  of  freedom  (r=9)  can  be  derived  directly  using  the  modified 
Hooker-Margulies  Equations  (6)  extended  by  Widhalm  and  Conway  (1)  to  take  into 
account  the  translation  of  Joint  1  relative  to  Body  0. 
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where 


^oo-ZL^'‘»  •  °  dyadic 


A  U 


°o*  =  ZZ^‘.-^‘«' .  a  vector 


K  u 


A  u 


=  scaiar 


A  « 


and 


6,^  =  1,  if  gi  belongs  to  a  joint  anywhere  on  the  chain 


of  bodies  connecting  (U  with  the  reference  body 


0,  otherwise  (e.g.  if  m  =  0) 


£r  =  25oi  +  2c0qX5q 


=  -Z?xZ?x]+  -25,^Z5,^ 


=  ^x^l  -Z5;,x^Xm] 


^X  =  -  z 


^  Km  ~  ^  K*  T-km 


where  is  the  vector  from  the  mass  center  of  a  body  a,  to  the  joint  on  body  A  leading 

to  a  body  n  ,  m^  and  <Px  are  the  mass  and  the  inertia  dyadic  of  body  a  ,  and  n  is  the  total 
mass  of  the  five-body  system.  The  actual  values  for  the  masses,  lengths,  and  inertia 

properties  used  in  the  model  are  specified  in  Chapter  IV.  Fi  is  determined  from 


K-Y.'^KM'Y.^KMykQK 
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and  Fi  is  the  vector 


^*k 

+  myp'^(l  -3pp)'  Z5^J 

In  Eq  (5). 

tOo  =  oJoiOi  +  a>02®2  <^0363 

and  the  y,  terms  again  refer  to  the  rotation  about  the  g,  axes  of  rotation.  Note  that 

ZqI  ~  ^02  “  ^03  “  ^04 

since  all  of  these  terms  represent  the  vector  from  the  center  of  mass  of  the  reference  body 
to  Joint  1,  the  translating  joint.  Note  also,  superscript  R  in  Eq  (15)  refers  to  the  fact 
that  the  time  derivatives  are  taken  relative  to  the  rotating  ei,  02,  63  reference  frame. 

Therefore,  the  DJ, ,  DJ,  terms  refer  to  the  translational  velocity  and  acceleration  of  Joint  1 
relative  to  the  surface  of  the  OMV.  The  spring  and  damper  torques  TSD  refer  to  the  six 
internal  motor  control  torques  TGI,  TG2,  TG3,  TG4,  TG5,  and  TG6.  The  external 
torques  T  refer  to  the  three  external  thruster  control  torques,  Tl,  T2,  T3  ,  located  on  the 
reference  body.  All  external  forces  are  assumed  zero.  In  addition,  all  gravity  gradient 

terms  are  assumed  negligible  (all  terms  containing  the  planetocentric  position  vector  p  ). 
The  joint  translational  motion  and  the  internal  and  external  control  torques  need  to  be 
specified  in  the  equations  of  motion,  which  is  the  topic  of  next  chapter. 

As  mentioned  earlier,  the  constraint  forces  and  torques  at  the  joints  have  been 
eliminated  from  the  equations  of  motion.  These  terms  can  easily  be  recovered  as 
described  in  references  (6)  and  (7).  The  equation  for  the  constraint  torques  is  as  follows, 
where  the  subscript  i  refers  to  the  particular  joint;  Joint  1,  Joint  2,  or  Joint  3: 


c 


-  ^'ij 

and  the  constraint  force  at  the  translating  joint  can  be  found  with 

^01  =  rrioTo 

where  to  is  the  position  vector  from  the  system  mass  center  to  the  center  of  mass  of  the 
reference  body.  Body  0. 
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There  are  two  types  of  control  schemes  required  in  this  thesis.  The  first  involves 
maintaining  dynamic  balancing  and  stability  throughout  the  deployment  of  the  grappling 
arm  and  the  subsequent  spin-up  of  the  grappling  device  prior  to  docking.  The  second 
control  problem  involves  the  detumbling  and  despinning  of  the  OMV/Target  docked  system 
to  a  spin-stabilized  state  of  equilibrium.  The  first  problem  incorporates  a  rather  simple  set 
of  open-loop  control  laws.  The  second  control  problem  is  more  complex  and  involves 
applying  a  Liapunov  analysis  to  the  system  in  order  to  derive  the  feedback  control  laws. 

Before  specifying  either  the  open-loop  or  the  closed-loop  control  laws,  a  quick  analysis 
should  be  performed  in  order  to  express  the  equations  of  motion  in  a  form  more  convenient 
to  work  with.  The  following  analysis  closely  parallels  the  formulation  performed  by 
Widhalm  and  Conway  (9:659-661). 

Eq  (5)  can  first  be  written  in  the  form 

Ax  =  F*  (2*^) 

where  A  is  defined  as  the  9  X  9  matrix  on  the  left-hand  side  of  Eq  (5).  The  vector  F*  is 
defined  as  the  nine  element  vector  on  the  right-hand  side  of  Eq  (5),  and 

x  =  [x,  Xj  ^3  •■^4  '^5  ^6  ^7  ^8 

—  [cOq]  (^02  *^03  Vl  y2  ”^3  ”^4 


X 


h  -hY 


(25) 


where  Eq  (25)  is  derived  from  the  state  variables 


X  =  [X,  X2  X3  X4  X5  Xg  X7  Xg  X, 

=  [c^01  «^02  ^^03  >1  72  -y3  'y4  ^ 


(26) 


Due  to  the  convenient  form  of  the  modified  Hooker-Margulies  equations,  the  control 
vector,  u  ,  can  easily  be  extracted  from  the  right-hand  side  so  that  Eq  (24)  can  now  be 
written  as 

.Tv  =  F  +  u  (27) 
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where 

U  =  [U,  U2  “3  “S  U.6  >J-7  “9]^ 

=  [71  72  73  7C1  TG2  TG3  7C4  7C5  7G6]^  (28) 

Tl,  T2,  T3  are  the  three  external  thrusters  and  TGI,  TG2,  TG3,  TG4,  TG5,  TG6  are 
the  six  internal  motor  torques  corresponding  to  the  six  axes  of  rotation  gi,  g2,  gs,  g4,  gs,  ge. 

The  symmetric  matrix  A  is  always  invertible  for  physical  systems,  and  pre-multiplying 
Eq  (27  )  by  the  inverse  of  matrix  A  yields 

X  =  +  A''u  (29) 

This  system  of  equations  of  motion  can  then  be  augmented  with  the  kinematical  equations 


X,,  *  X, 


X  =  X- 


X  *  X. 


where 


Therefore  ,  xio  is  the  angle  between  Body  1  and  Body  0,  xn  is  the  angle  between  Body  2 
and  Body  1,  and  xn  is  the  angle  between  Body  3/4  and  Body  1,  all  in  the  ei  direction.  X13 
is  the  angle  of  rotation  about  the  grappling  arm  itself.  It  is  used  only  in  the  coordinate 
transformations  required  to  rotate  the  inertia  matrices  to  the  ei,  e2,  63  system.  The  angles 
y,  and  y,  are  the  angles  of  rotation  of  Body  3  and  Body  4,  respectively,  about  the  gs  /  ge 
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axis  (See  Fig.  2).  Since  both  Body  3  and  Body  4  have  cylindrical  inertia  properties,  these 
two  angles  are  not  needed  in  any  of  the  coordinate  transformations  and  are  not  included  in 
the  kinematical  equations. 

This  completes  the  set  of  attitude  equations  of  motion.  The  augmented  state  vector  is 
now  defined  as  xia*  and  contains  the  original  state  vector  x  plus  the  additional  elements 

XlO>  Xll.  5ti2  ,  Xi3  . 

Note  that  Eq  (5)  rewritten  as  Eq  (27)  still  contains  the  terms 

Z  7^  7^ 

i-oi  >  r-oi  '  ^01 


corresponding  to  Joint  1  position,  velocity,  and  acceleration  relative  to  the  reference  body. 
Since  this  translating  joint  motion  is  a  specified  function  of  time,  Eq  (27)  is  a 
nonautonomous  system.  It  will  later  be  shown  that  it  is  desirable  to  work  with  autonomous 
systems  when  applying  the  Liapunov  analysis.  Widhalm  and  Conway  (9:660)  suggested  that 
the  joint  motion  can  be  expressed  as  a  third  order  linear  system.  The  joint  position 
relative  to  the  reference  body  is 

^01  ~  f  1  “  ^ 

Note  that  w  and  c  are  both  constants,  c  is  the  height  from  the  center  of  mass  of  the  base 
of  the  OMV  to  its  top  surface,  w  is  the  width  from  the  original  Joint  1  position,  when  the 
grappling  arm  is  stored  horizontally,  to  the  es  axis,  w  is  also  half  the  length  of  the 
grappling  arm  which  is  the  length  from  Joint  1  to  Joint  2  and  also  the  length  from  Joint  2 
to  Joint  3.  Yi  is  the  joint  translation  in  the  positive  e2  direction  (i.e.:  towards  the  6$  axis). 
Since  w  and  c  are  both  constants: 


^01  ^  1  ®  2  “  ^  2®  2 

Iqi  =  y  1 02  =  y  2^7.  ~  y  3®2 


(33) 

(34) 


Expressing  the  Joint  1  translational  motion  as  a  third  order  linear  system  and  combining  it 
with  Eq  (29)  and  (30)  yields  the  following  set  of  equations  for  the  system: 

x  =  A''T+A''u. 

X'  1 0  =  X  4 
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(35) 


12 


^  13  “  S 

y  =  DY 


where 


D  = 


0 

0 


1 

0 


0 

1 


^31  ^32 


D 
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(36) 


The  constants  D31,  D32,  and  D33  can  be  specified  for  whatever  Joint  1  motion  and  response 
time  is  desired. 


Open-Loop  Control  Laws 


As  previously  mentioned,  the  purpose  of  the  open-loop  control  laws  is  to  deploy  the 
grappling  arm  assembly,  and  to  spin-up  the  grappling  device,  while  the  OMV  system 
remains  in  a  state  of  pure  constant  spin  about  the  es  axis.  In  other  words,  wqi  and  wqz 
should  remain  zero,  and  W03  should  remain  constant.  In  addition,  the  HX  and  HY 
components  of  angular  momentum  should  also  remain  zero,  and  HZ  should  remain  constant 
if  the  external  thrusters  are  not  applied  (conservation  of  angular  momentum).  In  order  to 
maintain  this  state  of  pure  spin,  several  relationships  need  to  be  obeyed. 

Deployment  Phase 

At  all  times  throughout  the  deployment  process.  Joint  2  must  remain  along  the  63  axis. 
Also,  Body  2  must  form  a  "scissors"  arrangement  with  Body  1  so  that  the  63  axis  remains 
the  axis  of  symmetry  for  the  five  body  system.  Since  the  Joint  1  translational  motion  can 
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be  arbitrarily  specified,  the  various  requirements  for  the  dynamical  attitude  motion  of  the 
system  need  to  be  defined  as  functions  of  this  Joint  1  motion.  Since  Joint  2  must  always  be 
on  the  es  axis  the  following  relationship  has  to  be  satisfied; 


-^01  ^2  =  -  (>"1  -  «>)e2  =  «/cos(y ,  )e2 


so  that 


y  I  =  cos 


w  I 


Taking  the  first  and  second  time  derivatives  yields 

It/' sin!  y  I  !  tt/sinf  y  1  , 

and 

y'sSinfyi  i-rjyicosfy,; 

Vi  = - —Ti - ^ - 

ws\n  {  y ,  j 

_  y'3Sin(y,)-r2yiCos(y,) 
y;sin^(y,) 

Also,  the  "scissors"  arrangement  between  Body  1  and  Body  2  requires  that 

y3  =  -2y, 

>3  =  -2y, 

and 

y3^-2y, 
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(38) 


(39) 


(40) 


(41) 

(42) 


(43) 


Finally,  since  the  rotation  about  the  §2,  g^,  gs,  ge  axes  should  remain  fixed; 


I 


72  =  0 

y4  =  o 


(44) 

ys  =  o 

y6  =  o 


Eqs  (40),  (43),  and  (44)  specify  the  required  motion  of  Bodies  1,  2,  3,  and  4  during  the 
deployment  phase  in  order  to  maintain  dynamic  balancing.  The  required  control  torques 
now  need  to  be  calculated  in  order  to  achieve  this  motion. 

For  the  case  of  not  using  any  external  torques,  Tl,  T2,  and  T3  can  be  set  equal  to  zero. 
Referring  back  to  Eq  (29),  both  the  A  matrix  and  the  F  vector  are  known.  Also,  the  terms 

.V4,  Xg,  Xi,  X7,  Xg.  X% 

have  all  been  specified  in  Eq  (40),  (43),  and  (44).  The  only  unknowns  on  either  side  of 
these  six  equations  are  the  terms  contained  in  the  control  vector  u  .  Therefore  since  Tl, 

T2,  and  T3  are  zero,  the  only  unknown  terms  are  TGI,  TG2,  TG3,  TG4,  TG5,  and  TG6. 
There  are  now  six  equations  and  six  unknowns  which  can  be  solved  simultaneously  to  yield 
the  six  required  internal  control  torques. 

If  external  torques  are  allowed  to  ensure  the  system  remains  in  a  state  of  pure  spin, 
then  the  additional  requirements  that 

cijoi  =  0 
OUo2“*^ 

^  03  ~ 


can  be  added,  and  the  system  then  solved  for  all  nine  control  torques.  Rather  than 
simultaneously  solving  nine  equations  for  nine  unknowns,  these  nine  control  torques  can  be 
computed  directly  with  a  simple  manipulation  of  Eq  (27). 

u  =  Ax  -  T  (46) 
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It  is  desired  that  the  deployment  phase  be  accomplished  without  the  external  thrusters,  if 
possible.  For  either  case,  once  the  control  vector  has  been  solved  for,  it  can  be  substituted 
back  into  the  original  system  of  equations  which  are  then  numerically  integrated  to  yield 
the  state  vector. 


Spin-up  Phase 


This  same  methodology  can  be  applied  for  spinning  up  the  grappling  device  after 
deployment  of  the  arm  and  just  prior  to  docking.  Body  3  is  spun-up  to  match  the  spin  rate 
of  the  target  satellite  while  Body  4  is  simultaneously  spun- up  at  an  equal  and  opposite 
rate.  Since  Joint  1  is  no  longer  translating  and  is  stationary 

y.  =  o 

y3  =  0  (47) 


and  again 

y,  =  0 

y,  =  0  (48) 


but  this  time 


yg  =  constant  =  -y^ 


(49) 


Nonlinear  Feedback  Contol-Detumble/Despin  Phase 

Once  the  OMV  and  the  satellite  are  docked,  the  next  control  problem  is  to  detumble 
and  despin  the  target  by  applying  the  appropriate  internal  motor  torques  while  firing 
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external  thrusters  to  control  the  absolute  motion  of  the  system.  Joint  1  is  driven  to  the 
center  of  the  OMV  corresponding  to  the  63  axis,  and  Bodies  1,  2,  3  (including  the  target) 
and  4  are  all  erected  to  the  vertical  position.  The  OMV  should  remain  in  a  state  of  pure 
constant  spin  throughout  the  process. 

In  order  to  solve  this  problem,  a  feedback  control  approach,  similar  to  that  taken  by 
Widhalm  and  Conway,  needs  to  be  incorporated.  The  complete  nine  element  control  vector, 
u  ,  is  a  nonlinear  function  of  the  augmented  system  state  variables  and  the  Joint  1  motion 
expressed  in  Eq  (35).  Liapunov’s  direct  method  can  be  applied  to  derive  a  control  law 
which  is  globally  asymptotically  stable  with  respect  to  the  spin-stabilized  state  of 
equilibrium.  Eq  (35)  is  repeated  here  for  convenience. 

X  =  +  A'^u 

X  =  X4 


X  ,T=  X- 


r  =  DV 


Recall  from  the  previous  discussion  that  the  Joint  1  translational  motion  Y  is  expressed 
separately  as  a  third  order  linear  system  so  that  the  equations  of  motion  form  an 
autonomous  system  more  suitable  for  Liapunov  analysis.  Note  that  D  is  a  negative  definite 
matrix  selected  to  obtain  the  desired  decay  of  Joint  I  to  the  final  position  at  [0  0  c]. 

A  lemma  presented  by  Vidyasagar  (8)  ,  which  is  valid  for  autonomous  systems,  can 
now  be  applied.  It  is  stated  as  follows:  "Let  V(  X13*  ,  Y  )  be  continuously  differentiable 
and  suppose  that  for  some  d  s  0  the  set 

=  ,  7  :V[x\-.  ,  7]<d\ 


is  bounded.  Suppose  that  V  is  bounded  below  over  the  set  Sa*  and  that 

71  x].  ,  7]<0 
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for  all  xi3*  and  Y  in  Sa*  .  Let  S  denote  the  subset  of  Sa*  defined  by 

■S'“[x,3  ,  F  5'^;F(x:,3.F)-o] 


and  let  M  be  the  largest  invariant  set  of  a  system  which  is  contained  in  S.  Then  whenever 
xi3*  and  Y(0)  are  members  of  Sa*  ,  the  solution  of  the  system  of  Eq  (35)  approaches  M  as  t 
approaches  infinity."  (8:157) 

Since  the  system  of  Eq  (35)  is  autonomous  ,  every  state  trajectory  is  an  invariant  set. 
The  task,  therefore  ,  is  to  find  the  candidate  Liapunov  function,  V  ,  in  order  to  "derive  a 
nonlinear  feedback  control  law  that  drives  the  five-body  system  to  the  spin-stabilized  state 
of  equilibrium."  (9:660) 

On  the  same  lines  as  Widhalm  and  Conway  (9:660),  one  candidate  Liapunov  function  is 


where  I  is  the  identity  matrix,  Kjo,  K.ii,  K.i2  .  and  K13  are  positive  constants,  and  R  is  a 
positive  definite  constant  matrix.  The  function  is  continuously  differentiable,  and  taking 
its  derivative  with  respect  to  time  yields 

+  /:,3(x,3-90-’)x,3+/:,3x,3x.3  +  F"RF  +  F^Rr  (51) 

Substituting  from  Eq  (35)  gives 

F  =  X "/!  /I  -  ‘  F  +  /I  - '  u',  -  .Jx  ,0  -  90°)x,  +  /C , , ( X , ,  -  1  80" ) X, 

*  K X  K  ,,x  ,,x^*  V^[D^  R*  RD]7  (52) 

However  ,  since  R  is  positive  definite  and  D  was  specified  to  be  negative  definite,  the 
Liapunov  matrix  equation  can  be  written  as 


29 


Ik. V.'ft  ■"  &  t 


~ 


D'^  R  +  RD  =  -Q 

where  Q  is  a  positive  definite  matrix.  Therefore,  Y[DTR  +  RD]Y  is  negative  definite.  To 
make  the  derivative  of  V  at  least  negative  semidefinite,  choose  (9:660): 

u  =  -T 

-ABx  (54) 


substituting  back  into  (51)  yields 


+  +  /?d]7 


(55) 


If  the  matrix  B  is  positive  definite,  then  the  derivative  of  V  is  negative  semidefinite  in  Xio 
,  xii,  xi2.  However,  if  B  is  diagonal  with  positive  elements  (except  B33  *  0),  then  the 
derivative  of  V  is  negative  semidefinite  in  X3  and  Xio,  Xn,  xjj.  (9:660).  Note  that  90°, 

180°,  and  90°  have  been  subtracted  from  Xio,  x^,  X12,  respectively,  in  this  analysis.  This  is 
just  a  simple  coordinate  rotation.  The  spin-stabilized  state  of  equilibrium  is  about  the  X3 
axis  (e3  axis).  This  state  occurs  when  xio  =  90°,  xn  =  180°,  and  X12  =  90°. ,  so  that  the 
coordinate  rotation  is  needed  to  express  these  three  values  in  relationship  to  the  state  of 
equilibrium.  Whenever  these  three  values  are  not  90°,  180°,  and  90°  ,  respectively  (i.e.: 
when  they  are  nonzero  relative  to  the  X3  axis),  then  a  nonzero  control  u  results.  The  control 
law  for  u  in  Eq  (54)  represents  the  desired  nonlinear  feedback  control  law  to  drive  the 
system  of  Eq  (35)  to  a  spin-stabilized  state  of  equilibrium. 

Substituting  Eq  (54)  for  u  into  the  system  of  Eq  (35)  results  in  the  following  linear 
system  : 


.v-[0.0,0.-/C,o(.v,o-90‘’),-A:,3.v,3,-/:,,(.v,,-  180“1.  -/^,3(.v,3-90°),0,0  ' 


-Bx 


^  10  "  4 


Xi,  =  Xg 
■^12“  *^7 
^  1 3  “  S 

y  =  DY 


(56) 


Note  again  that  the  four  constants,  Kio,  Kn,  K12,  and  K13  are  positive  as  are  the 
elements  of  the  diagonal  matrix  B.  These  values  determine  the  motion  and  response  times 
of  the  various  bodies  as  they  decay  to  the  spin-stabilized  state  of  equilibrium,  just  as  the 
elements  of  the  D  matrix  determine  the  decay  of  Joint  1.  However,  since  the  rotation  (  y,  ) 
about  the  grappling  arm  itself  is  locked  prior  to  docking,  X13  and  xs  are  both  initially  zero. 
Therefore,  the  constants  K13  and  B55  will  have  no  effect  on  the  control  law  and  can  be 
chosen  arbitrarily.  Note  also  that  if  wqi  and  W02  (xi  and  X2)  are  initially  zero,  any  choice 
of  Bii  and  B22  will  still  result  in  these  two  components  of  angular  velocity  remaining  zero. 
Also,  B33  should  be  chosen  to  be  zero  so  that  W03  remains  constant.  The  choice  of  the  other 
K,  B,  D  constants  should  be  made  such  that  the  motion  of  Bodies  1,  2,  3,  and  4  is  closely 
coordinated  with  the  motion  of  Joint  1  so  that  Joint  2  and  the  target  center  of  mass  remain 
as  close  to  the  e3  axis  as  is  practical.  This  will  help  reduce  the  required  amount  of  control 
to  maintain  the  absolute  motion  of  the  system. 
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IV.  Results 

In  this  chapter,  the  mass  and  inertia  properties  for  the  five  bodies  of  the  OMV  system 
and  for  the  target  are  presented,  and  the  actual  docking  configuration  is  established.  The 
initial  conditions  for  the  three  phases  of  the  simulation  are  specified,  and  the  results  from 
each  phase  are  presented.  A  total  of  nine  cases  are  analyzed;  three  cases  from  the 
deployment  phase,  one  case  from  the  spin-up  phase,  and  five  cases  from  the 
detumble/despin  phase.  The  deployment  phase  is  run  without  external  thrusters,  with 
external  thrusters,  and  finally  with  external  thrusters  but  without  deployment  of  the 
counter-mass.  The  spin-up  phase  is  only  run  without  external  thrusters.  The  detum¬ 
ble/despin  phase  is  run  with  variations  in  response  time,  target  coning  angle,  and 
deployment  arm  length. 

The  equations  of  motion  formulated  in  this  problem  are  rather  complex,  and  the 
corresponding  computer  code  quite  lengthy.  To  validate  the  equations  and  computer  code, 
a  routine  was  added  to  compute  both  the  total  system  angular  momentum  and  the  total 
system  kinetic  energy  during  each  time  step  in  the  numerical  integration.  When  no 
external  torques  were  applied  (i.e.  Tl,  T2,  T3  =  0)  both  these  properties  were  successfully 
conserved,  thereby  verifying  the  validity  of  the  equations  of  motion  for  this  model. 

Both  the  target  and  the  OMV  mass  and  inertia  properties  are  presented  in  Table  I, 
and  the  necessary  dimensions  are  listed  in  Table  II.  The  target  properties  are  the  same  as 
those  used  in  the  Widhalm  and  Conway  study  (9).  The  base  of  the  OMV,  Body  0,  also  has 
the  same  properties  as  the  the  OMV  in  the  Widhalm  and  Conway  model.  These  values 
were  chosen  as  such  so  that  a  reasonable  comparison  could  be  made  between  the  results 
from  the  2-body  model  and  the  5-body  model.  An  estimate  was  made  for  the  mass  and 
inertia  properties  of  Body  1,  the  grappling  arm,  and  for  Body  3,  the  grappling  device. 

The  grappling  arm  was  chosen  to  have  a  total  length  of  3.5  meters,  with  two  equal  halfs 
of  1.75  meters  divided  by  Joint  2.  Body  1  includes  the  mass  of  the  translating  Joint  1 
assembly  so  its  center  of  mass  is  somewhat  to  the  left  of  Joint  2.  Body  4  has  exactly  the 
same  properties  of  Body  3,  and  together  they  form  two  equal  halfs  of  a  sphere.  Body  2’s 
required  inertia  properties  were  computed  using  Eqs  (2),  (3),  and  (4)  in  Chapter  2,  and  its 
corresponding  mass  was  estimated  given  these  requirements.  Note  that  the  mass  and 
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inertia  properties  of  Body  2  are  approximately  six  time  higher  than  those  of  Body  1.  This 
adds  considerable  extra  mass  to  the  system.  However,  the  total  mass  of  Body  1,  2,  3,  and 
4  is  390  kg  which  is  less  than  a  9%  addition  to  the  original  OMV. 


TABLE  I 

Nominal  Mass  Properties 

BODY 

MASS  (Kg) 

Ixx  (Kg-m) 

lyy  (Kg-m) 

I**(Kg-m) 

Target 

1000 

1000 

1000 

1100 

Body  0 

(Base  of  OMV) 

4500 

6400 

6400 

11800 

Body  1 

(Grappling  Arm) 

50 

50 

10 

50 

Body  2 

(Counter-mass) 

300 

311.125 

50.625 

311.125 

Body  3 

(Grappling  Device) 

20 

20 

20 

20.3125 

Body  4 

(Counter-mass) 

20 

20 

20 

20.3125 

TABLE  II 

Nominal  Dimensions  (Constants) 


L301  L31 


L401 


L30,  L31 


One  Half  the  Length  of  the  Grappling 
Arm  (i.e.  from  Joint  1  to  Joint  2  or 
from  Joint  2  to  Joint  3) 


Vertical  Length  from  Center  of  Mass  of 
the  Reference  Body  to  Joint  1 


Length  along  Grappling  Arm  from 
Center  of  Mass  of  Body  1  to  Joint  1 


Length  along  Grappling  Arm  from 
Center  of  Mass  of  Body  1  to  Joint  2 


Length  along  Grappling  Arm  from 
Center  of  Mass  of  Body  1  to  Joint  3 


(Pre- Docked)  Length  from  Center  of 
Mass  of  Body  3  to  Joint  3 


Length  from  Center  of  Mass  of  Body  4 
to  Joint  3 


(Post-Docked)  Length  from  Center  of 
Mass  of  Body  3  (Target  Satellite)  to 
Joint  3 


0.62  m 


0.35  m 


1.40  m 


3.15  m 


0.125  m 


0.125  m 


Docking  Configuration 

As  previously  mentioned,  the  target  used  in  this  problem  has  the  same  characteristics 
as  in  the  Widhalm  and  Conway  problem  (9).  These  characteristics,  along  with  the  target’s 
length  from  its  center  of  mass  to  Joint  3  after  docking,  determine  the  initial  positioning 

and  spin  rate  of  the  OMV  prior  to  deployment  of  the  grappling  arm,  and  the  necessary 

\ 

position  of  the  grappling  arm  assembly  after  deployment.  The  target’s  state  also 
determines  the  required  pre-docking  spin-up  of  the  grappling  device.  Nominal  motion  of 
the  target  is  listed  in  Table  III  and  its  mass  and  inertia  properties  have  already  been 
presented  in  Table  I. 
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The  initial  spin  rate  of  the  OMV  about  its  esaxis  of  symmetry  must  be  0.102  rad/sec 
in  order  to  match  the  precession  rate  of  the  target.  The  spin  rate  of  Body  3  must  be 
0.009  rad/sec  to  match  the  target’s  spin  rate.  The  docking  configuration  is  outlined  in 
Fig.  5.  It  can  easily  be  seen  that  with  the  given  lengths  and  coning  angles,  only  one 

configuration  will  work.  The  angles,  Vi  and  ,  required  for  docking  have  to  be 
determined.  Since  the  lengths  from  Joint  2  to  Joint  3  and  from  the  target’s  center  of  mass 
to  Joint  3  are  both  1.75  meter  ,  the  triangle  from  Joint  2  to  Joint  3  to  the  target’s  center 
of  mass  to  Joint  2  is  isosceles.  Therefore,  the  angle  formed  between  the  grappling  arm 

and  the  axis  is  equal  to  the  coning  angle,  20°  (0.3490  rads),  so  that  y,  must  be  70° 

(1.2217  rads)  in  the  docked  configuration.  Furthermore,  the  angle  >4  ,  the  rotation  of 
Body  3/Body  4  with  respect  to  the  grappling  arm,  must  be  50°  (0.8727  rads).  The 
distance  of  Joint  1  from  the  ea  axis  in  the  docked  configuration  determines  the  necessary 
amount  of  joint  translation.  This  distance  from  Joint  1  to  the  ea  axis  is  a  function  of  the 
target  coning  angle  and  is  equal  to  the  horizontal  displacement  ,  b  ,  of  Joint  3  from  the  ea 
axis.  This  distance  is  calculated  to  be  0.5985  meters. 

Deployment  Phase 

The  OMV  is  initially  pre-positioned  relative  to  the  target,  and  spun-up  to  match  the 
target’s  precession  rate.  Initial  conditions  for  the  deployment  phase  are  shown  in  Table 
IV.  Bodies  1,  2,  3,  and  4  are  all  stored  flat  on  the  top  surface  of  the  OMV’s  base.  Note, 
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however,  in  the  deployment  open-loop  control  law  Eq  (40)  there  is  a  singularity  at  y  i  = 

0°.  Body  1  is  therefore  stored  slightly  offset  from  the  horizontal,  at  an  initial  angle  of  1° 
(0.01745  rads),  and  Body  2  is  stored  at  an  angle  of  2°  (0.03490  rads).  This  is  a  simple 
work  around  to  this  singularity,  since  it  does  not  lend  itself  to  an  easy  removal  using 
L’Hospital’s  rule.  This  slight  offset  does  not  in  anyway  detract  from  the  analysis  of  the 
results.  250  seconds  are  allocated  for  the  total  deployment  event.  Joint  motion  is  specified 
as  a  third  order  linear  system,  with  three  equal  eigenvalues  of  -0.06,  giving  a  joint  motion 
of 

r.  =  -  0.000216fy,  +0.599) -0.0108^2-0. 18^3 
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TABLE  IV 

Deployment  Phase  Initial  Conditions 


0.0 

rad/sec 

0.0 

rad/sec 

0.102 

rad/sec 

0.01745  radians 

0.0 

radians 

0.03490  radians 

'  0.0 

radians 

0.0 

rad/sec 

0.0 

rad/sec 

0.0 

rad/sec 

0.0 

rad/sec 

0.0 

rad/sec 

■  ■ 

0.0 

rad/sec 

-1.75 

meters 

0.0 

m/sec 

0.0 

m/sec2 
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CASE  1  -  DEPLOYMENT  WITHOUT  EXTERNAL  THRUSTERS:  In  this  particular 
case,  the  deployment  of  the  arm  is  achieved  with  ihe  external  thrusters  turned  off.  It  can 
be  seen  from  Figs.  6,  7,  and  8  that  Joint  1  moves  smoothly  from  a  position  of  rest  -1.75 
meters  from  the  63  a.xis  and  decays  down  to  -0.5985  meters,  well  within  the  allocated  250 

seconds.  Looking  at  Vi  and  y,  in  Figs.  9  and  10,  the  grappling  arm  is  smoothly  raised 
from  loto  70°.  Body  2  is  successfully  deployed  (See  Figs.  11  and  12  )  in  a  manner  such 

that  Ys  is  always  twice  as  much  as  Vi  so  that  Body  2  forms  a  scissors  arrangement  with 
Body  1.  It  can  be  seen  from  Fig.  13  that  the  internal  motor  torque  TGI  responsible  for 
erecting  the  grappling  arm  assembly  experiences  a  sharp  start-up  transient  which  dies  out 
after  approximately  10  seconds.  The  peak  required  start-up  torque  is  less  than  7  N-M.  A 
sharp  transient  is  also  present  in  the  constraint  torques  at  Joint  1  and  Joint  2  (See  Figs.  19, 
20,  and  21  ).  The  rest  of  the  internal  motor  control  torques  (TG2,  TG3,  TG4,  TG5,  TG6) 
are  very  small  (Figs.  14-18  ),  since  they  all  play  minors  roles,  with  the  exception  of  TG3 
which  controls  the  motion  of  Body  2.  The  sharp  transient  control  torque  at  TGI  is  due  to 
the  fact  that  a  large  moment  is  initially  required  to  raise  the  arm  from  its  stored 
horizontal  position,  and  get  it  moving. 

The  most  important  goal  in  the  deployment  phase  is  to  maintain  dynamic  balance  and 
stability  throughout  the  process  without  having  to  fire  the  external  thrusters.  In  other 
words,  there  should  only  be  pure  spin  about  the  ea  axis,  and  wqi  and  wo2  should  remain 
zero.  Also,  the  cross  products  of  inertia  should  all  be  cancelled  out  by  the  deployment  of 
Body  2,  so  that  angular  momentum  exists  only  in  the  63  direction,  lined  up  with  wqs  . 

This  is  successfully  achieved.  Both  wgi  and  W02  remain  zero  throughout  the  process.  The 
HX  and  HY  components  (See  Figs.  23  and  24  )  are  negligible  when  compared  with  the  HZ 
component  in  Fig.  25,  and  can  be  attributed  to  numerical  round-off  errors.  Figs.  25  and 
26  show  that  all  the  angular  momentum  remains  completely  in  the  ea  direction.  Fig.  26 
also  shows  that  angular  momentum  is  conserved  throughout  the  process  since  no  external 
torques  are  applied,  thus  proving  that  the  equations  of  motion  for  the  system  are  valid. 

There  is  one  important  thing  that  does  not  work  as  expected,  however.  The  initial 
spin  rate  of  wga  =  0.102  rad/sec  does  not  remain  constant,  as  required,  but  instead 
increases  to  almost  0.106  rad/sec.  It  is  important  to  note  again  though  that  the  wqi  and 
W02  components  remain  zero  and  so  the  system  is  still  in  a  state  of  pure  ,  though  not 


39 


constant,  spin  about  the  es  axis.  The  reason  that  the  spin  rate  changes  is  due  to  a  simple 
conservation  principle  that  was  overlooked  in  the  initial  formulation  stages  of  the 
problem.  The  angular  momentum  remains  constant  about  the  es  axis  throughout  the 
deployment  process.  As  the  translating  joint  moves  inward  and  the  arm  is  erected,  the 
individual  angular  momentum  components  of  Body  1,  Body  2,  Body  3,  and  Body  4  all 
change,  but  the  total  s^^stem  angular  momentum  must  remain  the  same.  Therefore,  the 
angular  momentum  of  Body  0  must  somehow  change,  and  since  it  is  in  a  state  of  pure 
spin,  its  spin  rate  must  change.  This  is  not  a  problem  for  stability  but  it  is  a  serious 
docking  problem. 

The  spin  rate  of  the  OMV  must  match  the  precession  rate  of  the  target.  There  are 
two  possible  solutions  to  this  problem,  short  of  a  system  configuration  change.  The  first 
solution  merely  involves  applying  the  external  thrusters  during  the  deployment  process  to 
keep  the  spin  rate  of  the  OMV  fixed.  The  other  alternative  is  to  start  with  a  different 
initial  spin  rate  and  not  use  the  external  thrusters.  Since  the  desired  OMV  spin  rate  at  the 
end  of  the  deployment  process  is  known  (0.102  rad/sec)  and  angular  momentum  is 
constant,  the  required  spin  rate  at  the  beginning  of  the  deployment  can  easily  be 
computed.  This  was  done,  and  it  was  found  that  an  initial  OMV  spin  rate  of  0.0982 
rad/sec  will  yield  a  final  spin  rate  of  0.102  at  the  end  of  the  deployment  process  (i.e. 

when  yi  =  70°  and  Vs  =  140°)  .  This  lower  initial  spin-rate  saves  maneuvering  fuel  in  two 
ways.  First  of  all,  the  OMV  has  to  be  initially  spun-up  that  much  less,  and  second  ,  the 
external  thrusters  don’t  have  to  be  turned  on  at  all  during  the  deployment  process.  The 
main  problem  with  this  approach  is  with  alignment.  The  grappling  arm  assembly  is 
initially  lined  up  "underneath"  its  intended  docking  point  on  the  target.  If  the  OMV  does 
not  have  the  same  constant  spin  rate  throughout  the  deployment  process  to  match  the 
constant  precession  rate  of  the  target,  the  grappling  device  will  not  arrive  at  the  intended 
point  on  the  target.  If  the  OMV  is  spun- up  at  the  lower  initial  rate,  an  additional 
algorithm  is  required  to  ensure  that  the  grappling  device  arrives  at  the  intended  docking 
point.  This  adds  more  complexity  to  the  pre-positioning  program  and  runs  the  added  risk 
of  interference  with  appendages  on  the  target.  Therefore,  it  is  probably  more  desirable  to 
use  the  external  thrusters  to  keep  the  OMV’s  spin  rate  constant  about  the  es  axis  during 
the  deployment  process,  and  it  is  this  technique  that  is  incorporated  in  Case  2,  hereafter 
referred  to  as  the  Nominal  Deployment  Case. 
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One  last  noteworthy  piece  of  data  is  gleaned  from  this  case.  During  the  deployment 
process,  the  total  system  mass  center  of  the  OMV  assembly  remains  fixed  in  inertial  space 
since  no  external  torques  or  forces  are  applied.  However,  as  the  grappling  arm  assembly 
is  deployed,  the  system  mass  center  moves  "upward"  relativ..  the  base  reference  body. 

Since  the  system  mass  center  is  fixed  in  inertial  space,  the  base  reference  body.  Body  0, 
moves  "downward,"  away  from  the  target.  This  has  to  be  compensated  for  by  a 
translational  external  thruster  control  system  on  board  the  OMV  in  order  to  maintain  a 
constant  stand-off  distance.  This  amount  of  vertical  translation  that  needs  to  be 
compensated  is  computed  to  be  0.18  meters  as  shown  in  Fig.  27. 

CASE  2  -  NOMINAL  DEPLOYMENT:  The  second  deployment  run  is  made  using  exactly 
the  same  initial  conditions  as  Case  1  (Table  IV)  except  this  time  the  external  thrusters  are 
turned  on  in  order  to  maintain  the  OMV  at  a  constant  initial  spin  rate  of  0.102  rad/sec. 
Fig.  28  shows  that  the  control  scheme  is  successful  in  achieving  this  objective.  The 
internal  motor  torques  are  exactly  the  same  as  in  Case  1,  as  expected  since  their  task 
remains  unchanged.  The  required  external  thrusters  are  shown  in  Figs.  29,  30,  and  31  . 
Note  that  T1  and  T2  about  the  ei  and  e2  axes  are  both  negligible  and  can  be  attributed  to 
numerical  round-off  error.  The  only  external  control  torque  required  was  T3  about  the  03 
axis,  and  it  has  a  peak  magnitude  of  less  than  1  N-M.  Therefore,  the  amount  of  external 
control  required  to  maintain  a  constant  spin-rate  is  quite  reasonable 

CASE  3  -  DEPLOYMENT  WITHOUT  BODY  2  COMPENSATION:  In  both  Case  1  and 
Case  2,  Body  2  is  deployed  in  a  scissors  fashion  relative  to  Body  1.  This  successfully 
accomplishes  its  intended  purpose  of  cancelling  out  the  cross-product  of  inertia  terms 
thereby  maintaining  stability  without  the  use  of  the  external  thrusters.  However,  Body  2 
adds  a  considerable  amount  of  mass  to  the  system.  This  extra  mass  has  to  be  ferried  up 
into  orbit  along  with  the  OMV  and  also  has  to  be  transported  each  time  the  vehicle  is 
boosted  into  higher  capture  orbits.  This  could  amount  to  a  significant  fuel  cost.  There 
might  also  be  the  problem  of  Body  2  interfering  with  target  appendages  such  as  solar 
panels  and  antennas.  One  alternative  that  needs  to  be  examined  is  at  what  maneuvering 
fuel  cost  can  the  deployment  phase  be  achieved  without  deploying  Body  2.  In  this  case, 

Yo  is  fi.xed  ai  throughout  the  process  so  that  Body  i  and  Body  2  remain  coincident. 
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while  the  external  thrusters  are  turned  on  to  maintain  a  constant  spin  about  the  es  axis 
only.  The  resulting  control  magnitudes  are  shown  in  Figs  32,  33,  and  34  .  Most  of  the 
control  is  with  T1  about  the  ei  axis  in  keeping  the  system  balanced.  The  T2  and  T3 
thrusters  are  used  to  a  lesser  degree.  Both  their  magnitudes  are  less  than  1  N-M,  and 
neither  are  required  beyond  200  seconds  where  the  deployment  process  is  essentially 
achieved.  Notice  that  the  T1  thruster  is  still  turned  on  at  a  constant  l  N-M  at  the  end 
of  the  250  seconds  process.  This  appears  to  be  the  major  advantage  of  fered  by  Body  2. 
The  external  thruster  control  without  using  Body  2  is  not  unreasonable  during 
deployment,  but  must  be  kept  on  continuously  in  the  pre-docked  configuration.  By 
deploying  Body  2,  the  system  remains  in  a  state  of  pure  (although  not  constant)  spin 
without  the  use  of  external  thrusters.  The  advantages  offered  by  Body  2  versus  its  added 
mass  is  a  decision  that  would  need  to  be  weighed  by  the  mission  planner.  This  thesis 
merely  shows  that  both  options  are  feasible  from  a  controls  point  of  view. 

Spin- Up  Phase 

CASE  4  -  NOMINAL  SPIN-UP:  Once  deployment  is  completed.  Body  3  is  spun-up  in 
order  to  match  the  spin  rate  of  the  target  satellite  just  prior  to  docking.  Body  4,  the 
counter-mass,  is  spun-up  at  an  equal  rate  but  in  an  opposite  direction  in  order  to  cancel 
out  the  moments.  For  this  phase  of  the  simulation  Joint  1  is  held  fixed  as  are  the  other 

bodies  (wqs  =  0.102  rad/sec,  y  i  =  70« ,  and  -  140° ).  External  thrusters  are  shut 
completely  off.  Both  Body  3  and  Body  4  are  given  a  constant  acceleration  of  0.0009 
rad/sec2for  10  seconds,  and  then  the  OMV  system  is  left  in  its  pre-docked  pure  spin 
configuration  for  an  additional  40  seconds,  again  with  no  external  controls.  The  spin-up 
of  Body  3  and  Body  4  is  shown  in  Figs.  35  and  36  .  Both  bodies  successfully  spin-up  to  a 
constant  0.009  rad/sec.  Figs.  37  and  38  show  the  very  small  internal  motor  torques 
required  to  accomplish  this  task.  Figs.  39-42  show  that  dynamic  stability  is  maintained 
throughout  this  process  without  the  use  of  external  controls,  and  the  angular  momentum 
vector  remains  aligned  in  the  63  direction  along  with  the  constant  spin  rate. 


( 
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Detumble/Despin  Phase 


The  actual  mating  of  the  grappling  device  to  the  target  is  not  within  the  scope  of  this 
thesis,  and  occurs  as  a  discrete  event  in  which  Body  3  takes  on  the  mass,  length,  and 
inertia  properties  of  the  target.  Nonlinear  feedback  control  is  then  applied,  using  all  nine 
internal  and  external  controls,  to  drive  Joint  1  to  the  63  axis,  and  erect  all  bodies  in  the 

vertical  position  (yi  =  90° ,  >3  =  180°.,  >4  =  90° )  in  a  spin-stabilized  state  of  equilibrium. 
Also  in  this  process,  the  spin-rate  of  Body  3  (target)  and  Body  4  is  reduced  to  zero  in 
order  to  despin  the  target.  The  initial  conditions  for  the  detumble-despin  phase  are 
shown  in  Table  V. 


TABLE  V 

Detumble/Despin  Phase  Initial  Conditions 


0.0 

I 

rad/sec 

0.0 

rad/sec 

0.102 

rad/sec 

1.2217 

radians 

0.0 

1 

radians 

^  2.4434 

radians 

0.8727 

radians 

0.0 

rad/sec 

0.0 

rad/sec 

'  0.0 

( 

rad/sec 

'  0.0 

rad/sec 

!  0.009 

rad/sec 

-0.009 

rad/ sec 

-0.5985  meters 

0.0 

m/sec 

0.0 

m/sec2 

CASE  5  -  DETUMBLE/DESPIN  WITH  WIDHALM/CONWAY  GAIN  SETTINGS:  In  this 
case,  control  values  are  chosen  in  such  a  fashion  as  emulate  the  Widhalm  and  Conwav 
study  £9)  as  much  as  possible  for  comparison  purposes.  Feedback  is  applied  for  300 
seconds.  The  joint  translational  motion  is  specified  using  the  same  third  order  linear 
equation  as  Widhalm  and  Conway  (9),  corresponding  to  three  equal  eigenvalues  of  -0.04.  : 

/t  =  -0.000064K,  -0.0048ro-0.12rT  (58) 


The  constants  Kio,  Kn,  Kjj  ,  and  K13  as  well  as  the  nine  elements  of  the  positive 
diagonal  B  matrix  (See  Eqs  (54)  and  (56))  have  to  be  specified.  B33  is  set  equal  to  0.0 
since  it  is  desired  that  the  OMV’s  spin  rate  about  the  es  axis  remains  constant  throughout 
the  process.  B^  and  B22  are  set  equal  to  0.046  as  is  done  in  the  Widhalm  study.  Actually, 
since  wqi  and  W02  (corresponding  to  xi  and  X2  )  are  both  initially  zero  and  remain  that 
way  ,  the  constants  B^  and  B22  have  no  effect  on  the  problem  and  can  be  any  value. 

Note  also  that  the  rotation  about  the  g2  axis  (the  axis  coincident  with  the  grappling  arm) 
needs  to  be  initially  zero.  Rotation  about  this  axis  after  docking  and  prior  to  erection 
could  be  disastrous,  since  it  would  create  stability  problems  and  could  even  result  in  the 

target  crashing  into  the  base  of  the  OMV.  Therefore,  since  (X13  )  and  Vz  (xs  )  are 
initially  zero,  the  constants  K13  and  Bss  have  no  influence  on  the  problem.  K13  is 
arbitrarily  set  equal  to  a  value  chosen  for  one  of  the  other  K  constants,  and  B55  is  set 
equal  to  Bu  and  B22  • 

The  choice  of  B44  and  Kio  determines  the  second  order  linear  system  which 
represents  the  motion  of  y,  (xio)  from  70°  to  90°.  Similarly,  Bgg  and  determine  the 

motion  of  >3  (xn)  from  140°  to  180° ,  and  B77  and  K12  determine  the  motion  of  y4  (X12) 
from  50°  to  90° .  Widhalm  and  Conway  (9)  chose  to  specify  the  erection  of  the  target 
(which  is  similar  to  the  erection  of  Body  1  in  this  model)  as  a  critically  damped  system 
using  two  equal  eigenvalues  of  -0.035.  The  equivalent  motion  of  Body  1  in  this  problem 
needs  to  be  closely  coordinated  with  the  decay  of  the  translating  joint  in  order  that  Joint 
2  remain  as  close  to  the  es  axis  as  possible.  The  erection  of  Body  3  (the  target)  should  be 
timed  so  that  the  center  of  mass  of  the  target  also  remains  close  to  the  es  axis.  This 
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reduces  the  amount  of  external  torque  needed  to  maintain  a  dynamically  stable  system. 
Body  2  should  decay  in  a  fashion  similar  to  Body  1.  The  motion  of  all  three  angles  is 
chosen  to  be  represented  by  the  critically  damped  second  order  equation  consisting  of  two 
equal  eigenvalues  of  -0.035.  Therefore  B44,  Bge,  and  B77  are  all  chosen  to  be  0.07  while 
Kio,  Kii,  Ki2  and  are  all  chosen  to  be  0.001225. 

Finally,  the  decay  rate  for  the  spin  of  Body  3  (the  target)  is  determined  by  the  choice 
for  Bgg-  This  again  is  chosen  using  the  same  decay  rate  as  in  the  Widhalm  and  Conway 
model  (9)  ,  an  eigenvalue  of  -0.02,  making  Bgs  equal  to  0.02.  The  decay  of  Body  4  (m 
the  opposite  direction)  is  also  chosen  to  be  -0.02  so  B99  too  is  set  at  0.02. 

The  results  illustrate  that  the  detumbling  and  despinning  of  the  satellite  is  successfully 
accomplished  in  the  allotted  300  seconds,  leaving  the  system  in  the  desired  state  of  pure 
spin  shown  in  Fig.  4  .  All  motion  shown  in  Figs.  43-48  is  smooth  and  is  such  that  the 
target  center  of  mass  is  always  directly  above  Joint  2.  Body  2  has  exactly  twice  the  rate  as 
Body  1  at  all  times,  since  it  has  to  cover  twice  the  angular  displacement.  Fig.  49  shows 
the  successful  despin  of  Body  3  -  the  target,  while  Fig.  50  illustrates  the  corresponding 
despin  of  Body  4.  The  smooth  joint  translational  motion  is  shown  in  Figs.  51,  52,  and  53. 

The  most  important  information  from  this  run  is  the  required  external  and  internal 
control  torques  needed  in  the  detumble  and  despin  process.  The  largest  control  torques 
are  clearly  the  internal  motor  torque  TGI  and  the  external  thruster  Tl,  both  in  the  ei 
direction.  TI  starts  off  at  approximately  18  N-M  and  drops  down  to  almost  -9  N-M  in 
the  opposite  direction  before  decaying  towards  zero.  This  is  an  order  of  magnitude  higher 
than  the  external  torque  in  the  results  of  the  Widhalm  and  Conway  study  (9),  which  were 
bounded  from  +2  N-M  to  -3  N-M  . 

There  are  two  probable  reasons  why  the  required  TGI  and  Tl  and  also  TG4  control 
are  so  much  higher  than  in  the  Widhalm  and  Conway  model  (9).  First  of  all,  this  model 
has  the  additional  mass  and  inertia  of  the  grappling  arm  ,  the  counter-mass.,  the  grappling 
device,  and  its  counter-mass.  This  adds  390  kg  to  the  system,  which  is  a  39%  addition  to 
the  target  mass.  Second,  the  geometry  in  this  problem  is  considerably  different.  In  the 
Widhalm  and  Conway  model  (  Fig.  1),  the  target  was  joined  directly  at  the  translating 
joint.  In  this  model,  the  target  is  connected  at  Joint  3  which  is  linked  with  Joint  1  by  a 
3.5  meter  long  grappling  arm.  The  maximum  negative  control  torque  magnitude  for  both 
TI  and  TGI  occurs,  just  as  in  the  Widhalm  and  Conway  study  (9)  at  50  seconds,  and 


46 


corresponds  to  when  Joint  2  and  the  center  of  mass  of  the  target  both  have  the  equal 
maximum  displacement  from  the  63  axis  of  approximately  11.4  cm..  The  second  largest 
required  internal  motor  torque  is  seen  by  TG4  at  Joint  3,  which  is  responsible  for  erecting 

the  target  itself  relative  to  Body  1  (i.e.  y^).  The  amount  of  torque  required  to  despin  the 
target,  TG5,  is  quite  small  due  to  its  very  small  initial  spin  rate. 

All  internal  and  external  control  torques  decay  very  close  to  zero  within  300  seconds, 
implying  that  the  system  achieves  a  spin-stabilized  state  of  equilibrium.  The  magnitude 
of  the  constraint  torques  at  Joints  1,  2,  and  3  are  shown  in  Figs.  63,  64,  and  65  .  Notice 
that  the  magnitudes  of  these  torques  correspond  approximately  to  the  magnitude  of  the 
internal  motor  torque  at  the  particular  joints.  Also,  the  constraint  force  in  the  62 
direction  for  the  translating  joint  is  shown  in  Fig.  66  .  The  force  in  the  y-direction  is  the 
amount  of  force  required  to  "push”  the  translating  joint  towards  the  center  of  the  OMV. 

CASE  6  -  REDUCED  GRAPPLING  ARM  LENGTH;  The  next  case  examines  the  effect 
of  reducing  the  length  qf  thg  grappling  arm  to  see  if  a  shorter  arm  reduces  the  large  TGI 
and  T1  control  torques  seen  in  Case  5.  The  question  first  is  how  short  can  the  grappling 
arm  be,  and  still  accomplish  its  task.  Looking  at  Fig.  5,  it  is  apparent  that  with  the 

coning  angle  and  the  length  from  the  target  center  of  mass  to  Joint  3  fixed,  y<,  changes  as 
the  length  w  from  Joint  2  to  Joint  3  changes.  The  smaller  the  w,  the  less  the  Body  3  and 
Body  4  assembly  needs  to  be  rotated  in  order  to  dock  with  the  target.  As  a  design 

constraint,  it  is  decided  to  limit  the  minimum  angle  y*  at  0°  in  order  that  the  target  not 
run  the  risk  of  interference  with  the  grappling  arm  assembly  or  the  reference  body.  In 
other  words  when  docked  at  this  minimum  angle  of  rotation,  the  target  should  be 
perpendicular  to  the  grappling  arm  so  that  the  angle  between  the  length  from  the  target 
center  of  mass  to  Joint  3  and  the  length  from  Joint  2  to  Joint  3  is  90° .  The  length  w 
required  for  this  can  easily  be  computed  using  simple  geometry  ,  and  is  found  to  be 
0.63695  meters  so  that  the  total  length  of  the  grappling  arm  is  1.2739  meters.  This  is 

considerably  less  than  the  original  3.5  meters.  The  angle  Vi  is  also  found  to  be  equal  to 

the  coning  angle.  Therefore,  the  new  initial  conditions  are  yi  =  20° ,  y,  =  40° ,  and  y^  = 

0° .  Also,  the  joint  translation  and  the  the  control  values  for  K.  and  B  are  all  the  same  as 
those  used  in  case  5. 
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The  resulting  Tl,  TGI,  and  TG4  control  torques  are  shown  in  Figs.  67,  68,  and  69  . 
These  three  are  chosen  since  they  have  the  most  significant  magnitudes  in  Case  5. 
Shortening  the  length  of  the  grappling  arm  does  not  help  reduce  Tl  and  TGI.  In  fact, 
TG4  actually  increases  due  to  the  fact  that  it  rotates  the  target  from  0°  to  90°  in  the  same 
amount  of  time  it  originally  had  to  rotate  it  from  50°  to  90° . 

The  most  important  piece  of  information  gained  from  this  run  is  the  fact  that  TGI 
and  Tl  are  not  reduced  by  significantly  decreasing  the  length  of  the  grappling  arm  and 
instead  remain  almost  exactly  the  same.  There  is  one  thing  to  note  here  that  could  explain 
this  lack  of  change.  Even  though  the  length  of  the  arm  changes,  the  length  b  in  Fig.  5 
remains  exactly  the  same  .  Length  b  is  the  horizontal  displacement  of  the  end  of  the 
target  at  Joint  3  from  its  center  to  mass  on  the  63  axis,  and  is  a  function  of  the  target’s 
length  and  coning  angle  only.  The  length  b  is  also  the  length  that  Joint  1  has  to  translate 
inward.  It  appears  that  Tl  and  TGI  are  somehow  dependent  on  this  length  b. 

CASE  7  -  INCREASED  TARGET  CONING  ANGLE:  Since  altering  the  length  of  the 
grappling  arm  has  little  effect  on  Tl  and  TGI  ,  the  next  two  cases  study  the  effect  of 
altering  the  horizontal  length  b  mentioned  above  by  varying  the  initial  coning  angle  of  the 
target.  The  length  ,  precession  rate,  and  mass  properties  of  the  target  are  all  kept  the 
same,  but  the  coning  angle  is  increased  from  20°  to  30° .  This  decreases  the  initial  spin 
rate  of  the  target  to  0.008  rad/sec,  in  accordance  with  Eq  (1).  The  original  length  of  the 
grappling  arm  of  3.5  meters  is  used.  From  Fig  5  with  the  new  coning  angle  of  30°  it  can 

be  seen  that  the  new  initial  conditions  are  y,  =  60° ,  Vs  =  120°  ,  and  y4  =  30°  .  Joint 
translation  motion  and  decay  rates  are  again  the  same  as  in  Case  5.  This  leads  to  a  new 
length  b  of  0.875  meters  and  an  initial  Joint  1  position  of  -0.875  meters  from  the  the  ea 
axis. 

The  resulting  Tl,  TGI,  and  TG4  for  this  run  are  shown  in  Figs.  70,  71,  and  72  .  It  is 
apparent  that  increasing  the  cone  angle  thereby  increasing  the  required  amount  of  Joint  1 
translation  has  the  effect  of  increasing  the  required  Tl  and  TGI  control  torques. 

Case  8  -  DECREASED  TARGET  CONING  ANGLE:  This  next  case  examines  the  effect 
of  decreasing  the  coning  angle.  A  coning  angle  of  10°  is  used  which  results  in  initial 

conditions  of  y,  =  80°,  y3  =  160°,  and  y,  =  70° ,  and  an  initial  Joint  1  position  of  only 
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-0.30388  meters.  The  resulting  Tl,  TGI,  and  TG4  are  shown  in  Figs.  73,  74,  and  75  . 
Note  that  the  peak  control  magnitudes  are  reduced  to  less  than  10  N-M.  Therefore, 
reducing  the  coning  angle  significantly  reduces  the  required  control  torques,  TGI  and  Tl 
Case  7  and  8  illustrate  that  the  coning  angle,  which  determines  the  required  amount 
of  Joint  1  translation,  significantly  affects  the  principal  control  torques  Tl  and  TGI, 
while  Case  6  shows  that  altering  the  length  of  the  grappling  arm  has  very  little  effect. 

Case  9  -  NOMINAL  DETUMBLE/DESPIN:  The  next  and  final  case  incorporates  the 
original  grappling  arm  length  of  3.5  meters,  and  the  designated  target  coning  angle  of  20° 

.  It  uses  all  the  same  initial  conditions  as  in  Case  5  shown  in  Table  V,  The  big 
difference  though  is  th£  choice  q£  control  variables  KIQ.  K1 1.  K12.  and  the  elements  of 
the  diagonal  B  matrix.  Joint  translational  motion  is  again  kept  the  same  as  that  specified 
in  Eq  (58).  Bn,  B22,  B33,  K13,  B55,  Bgs  and  B99  are  also  chosen  to  remain  the  same  as  in 
Case  5.  The  values  that  are  manipulated  are  the  values  Kio,  B44,  Kn,  Bee,  K-i2>  B77  that 

specify  the  decay  functions  of  the  angles  Vi .  and  y*.  There  are  three  things  at 
stake  here  in  choosing  these  constants.  First  of  all,  these  constants  show  up  directly  in  the 
control  laws  (see  Eq  (54)),  so  that  higher  instructed  decay  rates  require  higher  control 
magnitudes.  Second,  decay  rates  determine  how  long  it  takes  the  complete  system  to  reach 

its  spin-stabilized  state  of  equilibrium  (yi  =  90°,  ys  =  180°,  and  y,  =  90°).  The  third  and 
final  effect  the  choice  of  these  control  constants  have  is  on  how  well  both  Joint  2  and  the 
center  of  mass  of  the  target  are  maintained  as  close  to  the  03  axis  as  possible.  It  is 

important  to  coordinate  the  decay  of  y,  and  y*  with  the  Joint  1  translation.  Large 
deviations  of  either  Joint  2  or  the  target  center  of  mass  call  for  large  control  magnitudes 
to  balance  the  system  and  keep  it  stable.  In  case  2  the  largest  negative  external  Tl  (and 
TGI)  control  torque  of  approximately  -8  N-M  occurs  50  seconds  into  the  detumble/despin 
phase,  when  both  Joint  2  and  the  target  center  of  mass  are  offset  by  almost  11.5  cm.  to 
the  left  of  the  e3  axis. 

Various  efforts  were  made  to  reduce  the  large  control  torques  Tl,  TGI,  and  to  a 
lesser  extent  TG4,  seen  in  Case  5  by  altering  these  response  values.  At  first,  the  decay 

rates  of  y,  and  yj  (Body  1  and  2)  were  kept  the  same  while  the  decay  rate  of  y^  (Body 
3/4)  was  altered.  Quicker  response  times  led  to  higher  initial  values  of  TGI  and  Tl. 


Slower  responses,  while  reducing  the  initial  control  magnitudes,  led  to  much  higher 
negative  torques  in  the  process,  due  to  the  fact  that  the  target’s  center  of  mass  became 
further  offset  from  the  es  axis.  It  appears  that  overall  lower  control  torques  are  produced 
when  Joint  2  and  and  the  target  center  of  mass  are  lined  up  vertically  as  in  Case  5. 
Therefore,  the  last  option  for  lowering  the  control  torques  is  to  change  the  original 
response  values  used  in  the  Widhalm  and  Conway  analysis.  Since  Joint  2  and  the  target 
cm  are  offset  to  the  left  of  the  ez  axis,  this  offset  can  be  reduced  by  slowing  the  decay 

rates  of  Vi .  and  while  keeping  the  joint  translation  the  same.  The  only  trade-off  is 
response  time.  The  final  values  chosen  are  equal  eigenvalues  of  -0.025  for  all  three 
functions.  This  leads  to  a  value  of  0.000625  for  Kio,  Kn,  and  K12,  and  a  value  of  0.05 
for  B44,  Bee,  and  B77.  Also,  due  to  the  slower  response  rate,  an  additional  50  seconds 
worth  of  feedback  is  applied,  allowing  a  total  of  350  rather  than  300  seconds  for  the 
detumble/despin  event. 

The  results  of  this  run  are  illustrated  in  Figs.  76  -  99  .  Body  1,  Body  2  and  Body  3/4 

(y  I .  Ys  •  yO  all  decay  smoothly  to  their  final  designated  positions  within  the  350 
seconds.  The  effects  of  changing  the  response  time  are  shown  in  Table  VI. 


TABLE  VI 

Effect  of  Decay  Rate  on  Final  Equilibrium  Position 


Case  5/  300  secs 


Case  9/  300  secs 


Case  9/  350  secs 


yi 

(degs) 

Vo 

(degs) 

Yi 

(degs) 

89.994 

179.987 

89.987 

89.906 

179.812 

89.812 

89.969 

179.938 

89.938 

It  can  be  seen  that  despite  the  slower  response  time  of  this  case,  the  bodies  still  come  very 
close  to  their  equilibrium  position  after  300  seconds.  The  additional  50  seconds  worth  of 
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feedback  control,  although  not  critical,  brings  them  even  closer  to  the  values  seen  in  Case 
5. 

Joint  1  translational  motion  and  the  despin  of  Bodies  3  and  4  are  the  same  as  in  Case 
5,  except  they  are  given  an  additional  50  seconds  to  damp  out  any  residue  motion.  The 
biggest  change  in  the  run  is  the  very  significant  improvement  in  the  TGI,  Tl,  and  TG4 
control  torques.  The  initial  Tl  and  TGI  torques  are  reduced  from  the  18  N-M  peak  value 
seen  in  case  5  to  a  much  lower  peak  value  9.31  N-M  .  Also,  the  negative  torque  is 
reduced  from  the  original  low  of  -8  N-M  at  the  50  seconds  point  to  a  much  smaller  value 
of  -1.63  N-M  at  the  43  seconds  point. 

The  reduction  of  the  initial  starting  torque  can  be  explained  by  the  slower  decay  rates 
fed  into  the  control  laws.  The  much  lower  negative  torque  at  43  seconds  can  be  explained 
by  better  coordination  of  these  decay  rates  with  the  Joint  1  translation  rate.  Joint  2  and 
the  target  cm  are  both  maintained  closer  to  the  ea  axis  of  symmetry.  In  case  5  they  are 
displaced  a  maximum  at  the  50  seconds  point  of  11.5  cm  to  the  left  of  the  es  axis.  In  this 
case  they  are  only  displaced  a  maximum  at  the  43  seconds  point  of  3.5  cm  to  the  right  of 
the  63  axis.  This  smaller  displacement  leads  to  much  smaller  required  control  torques,  Tl 
and  TGI.  The  TG4  torque  is  also  reduced  from  the  original  7  N-M  to  less  than  a  4  N-M 
initial  torque  with  almost  no  overshoot.  This  new  case  also  reduces  T2  and  smooths  out 
T3.  These  reduced  internal  and  external  motor  torques  provide  comparable  reductions  in 
the  constraint  torques 

The  results  of  this  case  clearly  show  a  significant  improvement  in  the  required 
external  and  internal  control  torques.  Reducing  the  response  rates  makes  a  big 
improvement,  while  only  having  to  extend  the  required  maneuvering  time  shortly. 
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Fig.  8.  Deployment  Phase  /  Case  1 
Joint  1  Relative  Acceleration 
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Fig.  10.  Dct  '  <vment  Phase  /  Case  1 
Booy  1  yi  History 
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Fig.  12.  Deployment  Phase  /  Case  1 
Body  2  Ya  History 
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Fig.  14.  Deployment  Phase  /  Case  1 
Internal  Motor  Torque  TG2 
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Fig.  16.  Deployment  Phase  /  Case  1 
Internal  Motor  Torque  TG4 
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Fig.  17.  Deployment  Phase  /  Case  1 
Internal  Motor  Torque  TG5 
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Fig.  19.  Deployment  Phase  /  Case  1 
Constraint  Torque  at  Joint  1 
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Fig.  25.  Deployment  Phase  /  Case  1 
Angular  Mome.:tum  Component  HZ 
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VERTICAL  TRANSLATION  OF  OMV 
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Fig.  27.  Deployment  Phase  /  Case  1 
Vertical  Translation  of  OMV 
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I  •  Fig.  30.  Deployment  Phase  /  Case  2 
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Deploymeiit  Phase  /  Case  3 
External  Torque  T1 
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Fig.  33.  Deployment  Phase  /  Case  3 
External  Torque  T2 
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Fig.  35.  Spin-up  Phase  /  Case  4 
Body  3  ys  History 
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NTERNAL  MOTOR  TORQUE  TG5 


Fig.  37.  Spin-up  Phase  /  Case  4 
Internal  Motor  Torque  TG5 
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Fig.  38.  Spin-up  Phase  /  Case  4 
Internal  Motor  Torque  TG6 
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Fig.  41.  Spin-up  Phase  /  Case  4 
Angular  Momentum  Component  HZ 
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Fig.  42.  Spin-up  Phase  /  Case  4 
Total  System  Angular  Momentum 
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Fig.  46.  Detumble/Despin  Phase  /  Case  5 
Body  2  y.  History 
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Fig.  49.  Detumble/Despin  Phase  /  Case  5 
Body  3  Ys  History 
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Fig.  51.  Detumble/Despin  Phase  /  Case  5 
Joint  I  Position 
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INTERNAL  MOTOR  TORQUE  TGI 


Fig.  57.  Detumble/Despin  Phase  /  Case  5 
Internal  Motor  Torque  TGI 
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Fig.  59.  Detumble/Despin  Phase  /  Case  5 
Internal  Motor  Torque  TG3 
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NTERNAL  MOTOR  TORQUE  TG4 


Fig.  60.  Detumble/Despin  Phase  /  Case  5 
Internal  Motor  Torque  TG4 
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Fig.  63.  Detumble/Despin  Phase  /  Case  5 
Constraint  Torque  at  Joint  1 
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Fig.  65.  Detumble/Despin  Phase  /  Case  5 
Constraint  Torque  at  Joint  3 
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Fig.  69.  Detumble/Despin  Phase  /  Case  6 
Internal  Motor  Torque  TG4 
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Fig.  70.  Detumble/Despin  Phase  /  Case  7 
External  Torque  T1 
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Fig.  71.  Detumble/Despin  Phase  /  Case  7 
Internal  Motor  Torque  TGI 
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Fig.  73.  Detumble/Despin  Phase  /  Case  8 
External  Torque  T1 
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Fig.  74.  Detumble/Despin  Phase  /  Case  8 
Internal  Motor  Torque  TGI 
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Fig.  75.  Detumble/Despin  Phase  /  Case  8 
Internal  Motor  Torque  TG4 
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Fig.  76.  Detumble/Despin  Phase  /  Case  9 
Body  1  Yi  History 
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Fig.  80.  Detumble/Despin  Phase  /  Case  9 
Body  3/4  History 
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Fig.  84.  Detumble/Despin  Phase  /  Case  9 
Joint  1  Position 
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Fig.  87.  Detumble/Despin  Phase  /  Case  9 
External  Torque  T1 
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Fig.  96.  Detumble/Despin  Phase  /  Case  9 
Constraint  Torque  at  Joint  1 
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Fig.  98.  Detumble/Despin  Phase  /  Case  9 
Constraint  Torque  at  Joint  3 
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V.  Conclusion 

The  equations  of  motion  were  derived  for  the  five  rigid  body  system  used  to  model 
capture  of  a  freely  spinning  and  processing  axisymmetric  target  satellite.  Control  laws 
were  formulated  for  each  of  the  intended  phases;  deployment,  spin-up,  and  detumble/de¬ 
spin.  A  Liapunov  analysis  was  performed  for  the  nonlinear  feedback  control  required  in 
the  detumble/despin  phase,  in  order  to  drive  the  system  to  a  spin-stabilized  state  of 
equilibrium.  A  computer  simulation  was  then  used  tc  examine  each  of  the  three  phases  of 
the  capture  process. 

The  deployment  of  the  grappling  arm  and  its  counter-mass  was  successfully 
accomplished  within  an  allotted  250  seconds.  A  sharp  transient  was  initially  observed 
with  the  TGI  motor  torque  at  Joint  1  due  to  the  large  moment  required  to  start  the  arm 
moving  from  its  stored  position.  However,  the  momentary  peak  was  less  than  7  N-M  and 
quickly  leveled  out.  It  was  shown  that  dynamic  balancing  and  stability  could  be 
maintained  without  the  use  of  any  external  torques,  but  that  the  initial  spin  rate  about  the 
63  axis  of  symmetry  did  not  remain  constant  due  to  the  conservation  of  angular 
momentum.  A  small  application  of  the  external  T3  thruster  about  the  63  axis  was  able  to 
alleviate  this  problem  and  keep  this  spin  rate  constant  as  required  for  capture. 

The  necessary  external  control  torques  to  achieve  deployment  without  utilizing  the 
counter-mass  were  also  shown,  in  order  to  illustrate  this  as  a  feasible  alternative. 

Once  the  grappling  arm  was  deployed,  the  grappling  device  was  spun-up  for  10 
seconds  and  its  counter-mass  was  spun-up  at  an  equal  and  opposite  rate.  An  additional  40 
seconds  were  allowed  to  elapse  to  observe  the  OMV  system  in  a  state  of  pure  spin  just 
prior  to  docking.  This  entire  50  second  process  was  achieved  without  any  external 
torques,  and  the  rysterfi  remained  in  its  state  of  pure  constant  spin.  The  spin-up  of  the 
grappling  device  and  the  counter-mass  required  very  small  internal  motor  torques. 

Docking  was  assumed  to  be  a  discrete  event  that  initiated  the  detumble/despin  phase. 
Nonlinear  feedback  control  was  immediately  applied  in  order  to  detumble  and  despin  the 
target.  With  the  final  selection  of  response  values,  feedback  was  applied  for  350  seconds. 
All  motion  of  the  bodies  was  smooth  as  they  travelled  to  the  spin-stabilized  state  of 
equilibrium.  The  response  values  were  chosen  so  that  the  movements  of  the  bodies  were 
coordinated  with  the  decay  of  the  translating  joint  to  ensure  that  the  composite  center  of 
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mass  stayed  close  to  the  es  axis.  The  most  significant  control  torques  observed  were  with 
the  TGI  internal  motor  at  Joint  1  and  the  external  T1  thruster,  both  about  the  es  axis. 
Each  initially  were  almost  10  N-M,  which  is  more  than  three  times  the  control  torques 
observed  in  the  Widhalm  and  Conway  two-body  model.  These  higher  torques  can  be 
explained  by  the  additional  mass  and  inertia  added  at  Joint  1  by  the  grappling  arm 
assembly  and  counter-mass,  and  by.  the  different  geometry  of  this  problem  with  the 
grappling  arm  acting  as  a  linkage  between  Joint  1  and  the  target. 

Additional  cases  were  examined  in  the  detumble/despin  process  for  design  purposes. 
Quicker  response  times  led  to  higher  initial  control  torques.  Reducing  the  length  of  the 
grappling  arm  had  no  significant  effect  on  control  magnitudes.  What  did  have  an  effect 
on  these  control  torques  was  the  change  in  coning  angles.  The  amount  of  required  control 
torques  clearly  increases  with  increasing  coning  angles.  Therefore,  control  requirements 
for  the  detumble/despin  process  are  direct  function  of  coning  angle.  It  was  also  shown 
that  torque  requirements  can  be  reduced  by  effective  control  of  the  the  target’s  center  of 
mass  position  relative  to  the  OMV’s  axis  of  symmetry.  These  are  all  important 
consideration  in  the  design  of  an  orbital  capture  vehicle. 

All  major  objectives  were  accomplished  in  this  thesis.  The  simulation  showed  that  a 
5-body  OMV  using  a  combination  of  internal  and  external  thrusters  can  successfully 
capture  and  detumble/despin  an  axisymmetric  target  satellite.  Counter-balance  masses  can 
be  effectively  used  in  lieu  of  external  thrusters  during  deployment  and  spin-up.  The 
five-body  system  definitely  adds  both  realism  and  flexibility  to  the  original  two-body 
model.  This  study  gives  a  good  indication  of  the  dynamics  and  controls  involved  in 
remote  orbital  capture  using  the  proposed  five-body  OMV  model,  and  has  hopefully 
paved  the  way  for  possible  future  studies  of  actual  hardware  design. 


There  are  several  possibilities  for  follow-on  work  with  the  OMV  model  developed  in 
this  thesis.  One  option  would  be  to  add  an  internal  momentum  wheel  about  the  63  axis  to 
keep  the  OMV’s  spin- rate  constant  during  deployment,  without  having  to  fire  the  T3 
thruster.  A  further  extension  of  this  idea  would  be  to  add  a  complete  set  of  three 
orthogonal  momentum  wheels  instead  of  external  thrusters  for  attitude  control,  as  was  done 
by  Fleming  (3)  on  the  Widhalm  and  Conway  two-body  model.  Another  possible  area  of 
follow-on  research  could  be  to  examine  a  wider  range  of  targets  and  capture  scenarios. 

This  particular  study  was  limited  to  axisymmetric  target  satellites  only.  The  capture  of  an 
asymmetric  target  would  present  additional  problems  and  requirements  that  would  have  to 
be  examined. 
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